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national advisory committee for aeronautics. 


TECHNICAL MEMORANDUM NO. 335. 

PRESSURE DISTRIBUTION ON JCUKOWSKI WINGS.* 

By Otto Blumenthal. 

In the winter semester of 1911-12, I described., in a lec- 
ture on the hydrodynamic bases of the problem of flight, the 
potential flow about a Joukowski wing.** In connection with 
this lecture, Earl Toepfer and Erich Trefftz computed the 
pressure distribution on several typical wings and plotted 
their results. I now publish these diagrams accompanied by a 
qualitative discussion of the pressure distribution, which 
sufficiently indicates the various possible phenomena. For a 
quicker survey, I have divided the article into two parts, 
the first part dealing with the more mathematical and hydrody- 
namic aspects and the second part, which is comprehensible in 
itself, taking up the real discussion from the practical stand- 
point . 


* From "Zeitschrift fur Flugtechnik und 
Hey 31, 1913. 

** See above magazine, Vol. I (1910), p. 


Motorluf tschif fahrt, " 
381. 
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I 

We obtain the entire number of all Joukowski wings of the 
length 2 1 with the trailing edge at the point x ~ - l, by 
laying, in a £ = t + i h plane through the point £ = l/2, 
the cluster of all the circles which contain the point £ = i/2, 
either inside or on their circumference, and plotting these 
circles by means of the formula 

Z = t. + IL (l) 

on the z . = x + i y plane- The circles, which contain the 
point' £ ~ i/2 on their circumference, thus become doubly in- 
tersected arcs and, in particular, the circle, which has the 
distance (-1/2, + i/2) for its diameter, becomes the recti- 
linear distance of the length 2 l . The circles which contain 
the point £ = l/ 2 inside, furnish the real Joukowski fig- 
ures. The point £ = - 1/ 2 passes every time into the sharp 
trailing edge- The individual Joukowski wings are character- 
ized by the following quantities (Fig. l). The center M of 
the circle K is connected with the point H, £ = - l/2, and 
the point of intersection of this connecting line with the T] 
axis is designated by 11*. The distance OM* on the *1 axis 
is equal to half the height of the arc produced by describing 
the circle about M* as its center and is therefore designated 
by f/2, as half the camber of the Joukowski wing, f being 
its first charact eristic dimension. We have chosen as the 
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second characteristic dimension, the radii difference Mil* = 6* 
This gives a measurement for the thickness of the Joukowski 


wing. 


We will now consider the determination of the velocity 
and pressure distribution which produce an air flow along the 
wing, in infinity, with the velocity V at an angle of n-j3 
with the positive x axis, 3 being the angle of attack of the 


wing. 


The absolute velocity q of this flow is calculated thus: 
If K (£,t]) is the absolute velocity of the air flow, of ve- 
locity V and angle of attack (3, around the circle K in 


the £ olane, then 


q (x,y) 




as i 

'd^l 


It is, however, 


ds ; 

d£l 


\ /(V - 4^; + I s Tl 2 , 

a j \ ** / 


a 2 = 


e 2 + n s 


and along the circle K 


k (I ,b) 


V~T 5 ~T 


+ 6 


3T(£ sinfB + T] cos 6) + c 


where 2 tt c is the circulation. This constant is determined 
according to Kutta, by the condition that the velocity at the 
trailing edge is finite and therefore, since dz/d£ there dis- 
appears, k must also disappear at the point H. Thus we ob- 
tain 
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K ( t , r , ) 


2 ft + T') sinP + r\ cospl. 


I / . -i \ j > 

Tv + tj sin ^ + ^ cos ! 3 ! (2) 


J ij + f ; i 


J v 


? 2 

fa l 2 - -L') + T 

v. 4 / 


a 2 = £ 2 + "T • J 

This rather involved expression is simplified by the in- 
troduction of a new v ariable, the angle w at the center of 
the circle K, measured from the radius MH. In this angle, 
tile coordinates t , rj and the quantities connected with them 
are expressed as follov/s: For abbreviation, v;e designate the 

tj i r2 

radius of the circle K with r = — ~ — =— + 6 and introduce 

the angle A by 


cos A - - 


f + r 


-, sin A = 


J I s + f 2 


The geometric significance of A and oj is obvious from Fig. 1. 
3y si mole calculations we now obtain 


1 - - \ + 2 r s in ^ s in + a'') , 

2 2 \ 2 / 


2 r sin “ cos + a'}, 

tO *■ i-u / 


— - 21 r sin — sin - ^ + A' + 4 r 2 sin 2 — 
2 \2 / 2 


l 2 , , 0) f c . w f /(JO .\" 

— + 4 r sm - I 6 sm - - - cos ( - + A ) 

t: I Cj fa Cj / 
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Formula 2 for q is simplified by the introduction of the 
angle co , to 



From this re next derive a few general results which, hold 
good for all the quantities f, 6 , (3 . 

a) On the trailing edge £ = — cos (a + 8). 

V 2r 

b) On top of the wing, there is always a portion along 

vjhich the velocity q > V, hence where there is a negative 
pressure. As proof of this, we will consider the center of 
the upper side, the point - —■ - a* At this point a > r, 

id 

as can easily be seen geometrically (Fig. l) or from formula 4. 

3tt 

However, if we put <*> = — - A in formula 5, it then becomes 

id 


q 

v 


CO£ 


(H - A _ 

V 4 2 H ' 


COS . T - 7? 


Ah 

3 / 




1 + 


H 


5 + 


Tf 

2/ 


COS 


/ tt 


cos 


-4-fh 

£— 3 J_ 

(TT _ Ah 
\4 2 j 


J 


7 T 

* A 


( 6 + 

V 3 / 


TT 


cos 


cos 


f TT _ 

13__ 

/ TT 


a _ 

. 1 lL 

A'-, 

2 / 
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c) The velocity is zero at the point u) = tt - 2a - 23, 
whicn is alv/a ys located on the lower side- At this point the 
streamline enters the wing. Further general conclusions (i-c., 
applying to all f, 6, [3 ) can hardly be drawn. We obtain 
considerably more accurate expressions in the especially int- 
eresting practical case where, in the vicinity of the leading 
edge, a pronounced velocity maximum and consequently a strong 
suction is produced- We will confine ourselves to this case 
in all tnat follows. Hereby we can, in formula 5, first of 
all disregard the slight fluctuation of the factor a 2 for 
small values of f and 6 and consider only the factor F, 
which must be alone decisive for the great changes in velocity. 
This factor, however, enables a simple explanation. 

For this puroose, we introduce the angle = jg- + A + 3 . 
The entering point of the streamline then lies at i| f = tt/2, 
where F disappears. In general, we have 


x 

f k 


a = 


(a 2 + sin 2 3) tan 2 t - 2 (ab - sin3 cosP) tani + 

+ (b 2 + cos 2 (3 ) , 

f (6 cos (A + 3 ) - ■£ sin(3 ^ , 


( 6 ) 


b = y- f 6 sin (A + p) + — cos(3^ . 


Consequently, F attains its maximum value at the angle ty 0 , 
which is given by the formula 


tan \|r 0 = ab sin . B cos(3 
a 2 + s in 2 |3 


(7) 
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and 


this value is 


max 


J ( 6 cos (A t- 


f Np 

} - o sin 6' + - — sin 2 (3 

nj </ “ 


5 cos A 


( 7 «) 


We now make the assumption, corresponding to the already 
announced purpose of our investigation, that F has a high 
maximum in relation to the value of cos (A + 3) on the trail- 
ing edge* We require, e.g. , that F rax shall equal or exceed 
J ?,. This is mathematically the most favorable. Formula 7', 
with the rid of a rough estimate, then gives 

Sin6 $ 0 COS (A + (?) [ 1 + — - p=~==r=r . 

2 v / I s + f* 

= 6 cos (A + 6) (1 + sin A) (3) 

With this insertion, the numerator of tan i|r Q is smaller than 

] j j — ~ 2 

- cos (A + p) | — — — — — co s j? - 6 s in ( A + p ) ( 1- sin A) 

r u 2 j 

For small f, 6 and 5, this value is always negative and there- 
fore the maximum value of F is assumed to be at a point lo- 
cated between the entering point and the trailing edge on the 

portion of the surface belonging to the upper side-* -On the 

* Generally the ooint is located on the upper side- It lies 
between the entering point and the leading edge, only when 6 
is very small in comparison with f. For 6=0, it lies on 
the leading edge. 
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other hand, it can be shown that the maximum is located not 
far from the entering point. In fact the greatly preponderat- 
ing member in the numerator of tan \|/ Q , on account of formula 
8, is cos0 sin (3 . The ease is not quite so simple with the 
denominator, which is 

a 2 + s irP 6 " 

= b 2 cos 2 (A + 0) - 6 f cos (A + 0)sin6 + — — - + - ■ sin 2 0- 

VI 3 - ^ 4 J 

If we introduce into the first member, on the right side of 

formula 8, the above limit for 6, the denominator is then 

smaller then 2 — — sin 2 p • Hence tan is either smaller 

1 

or at most only unessentially* greater than - i cotB, which 
shows that is either smaller or at most only slightly great- 
er than ■- + 2 0. The point go, at which F assumes its max- 
imum value, is located between the entering point of the stream- 
line and the upper side and, at most, only slightly farther 
than 4 0 from the entering point. 

Lastly, it may be remarked that in formula 6 for -q=r, both 

F ' 

the powers, tan 2 '!' and tanty, appear to be multiplied by 
small coefficients. Both these members therefore assume quite 
large values for large values of tan , i.e. , in the immedi- 
ate vicinity of the entering point and the leading edge. Hence 
F differs but little over the whole surface, with the excep- 
tion of the specified region, from the value cos (A + (3; on 

the trailing edge and th erefore only slightly f r om unity. 

* "Unessentially" means that the deviations are of the order of 
magnitude f/l , b/l . 
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The results of the computations in I arc as follows: A 

Joukowski wing is characterized by the three dimensions, name- 
ly, the length 21 , the camber f and the radii difference 
6, which is expressed 'in the thickness of the wing. To these 
is added the angle of attack jo . The points on the surface 
are most conveniently computed with the aid of a variable oo of 
the angle at the center of the circle in Fig. 1. The formulas 
for the coordinates will not be given here- In practice, a 
graphic process is employed which is explained in the accompa- 
nying note by E • Trefftz (pages 130-131 of this same volume 
of "Zeitschrift fftr Flugtechnik und hotorluftschiffahrt . 11 ) 

We require only the following data: w = 0 gives the trailing 

edge; oo = tt - 2A (tan A = ■— ) gives the leading edge and the 
intermediate values of go correspond to the lower surface of 
the wing. w =TT -2A-2P gives the entering point, i.e., 
the point where the air flow strikes the surface and hence 
where the velocity is zero. 

The ratio of the absolute velocity q of the air flow 
on the wing to the velocity V in infinity is given by 


q - ' jS tp 

v - ~T F 
r 2 


f + f 2 


is the radius of the circle in Fig. 1 and 
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a the distance OB in the same figure. Both factors, a and 
F, depend on u> . For small f and 5, the value of the factor 

< j 2 

— - differs but little from unitv. The orooerties of the fac- 

r I ' ‘ 

tor F, as obtained by the calculations of I, can be summarized 
as follows. 

On the trailing edge, F lias the value cos (A + (3) and 
decreases at the customary angles of attack (about 6°), from 
the leading edge to the entering point, where it becomes zero. 
For small f and 6, the decrease takes place very slowly 
throughout most of the lower side and first becomes rapid in 
the immediate vicinity of the entering point. 

From the entering point, F increases rapidly and attains 
near the leading edge, a maximum of the order of magnitude 
pi / 23. This is approximately also the maximum value of the 
velocity ratio q : V. For this maximum value, the ratio of 
the angle of attack to the thickness of the wing is therefore 
decisive, the camber having, in the first order, no effect on 
it. In constructive wing shapes, where 6 and P arc of the 
same order of magnitude, the velocity at the leading edge is 
accordingly not very great. This result is important because 
it explains the effect of rounding the leading edge* The re- 
sult is still more striking when we consider the radius of 


curvature p of the leading edge* It is, namely, with unes- 
sential omissions £. = 16 -4- ( 1 - 4 — \ The radius of curvature 

l f- \ IS 

therefore diminishes rapidly with decreasing 6, the rounding 
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off of the leading edge being very slight, and the maximum ve- 
locity remains within moderate bounds. The negative pressure 
on the leading edge, which, according to 3ernouilli’s equation, 
is proportional to q 2 /V 2 , is computed by the introduction of 
the radius of curvature in the first approximation, to 
4 p 2 I consider this simple formula worthy of attention. 

The course of F e.long the top of the wing can finally 
be characterized as follows: At some distance from the leading 

edge, F changes but slowly. If, therefore, the maximum value 
of F is much greater than unity, it falls abruptly at first 
and then gradually approaches the value at the trailing edge- 

Only in the vicinity of the leading edge does the factor 
F give us sufficiently accurate information concerning the 
course of the velocity q. Everywhere else we need to know 
th'e course of o . This can be easily found geometrically from 
Fig- 1. On the leading edge a lias the value of l/2. From 
the triangle HMO, it follows that a assumes its minimum value 


for the angle to 


min ’ 

sin to 


which is given by 


mu 


Jt 

2 


sin A V 5 2 cos 2 A + x sin" A 


Hence 


;m oo 


1 


mm 


/ + 4jl 


V f" f 

For the angl e 2 oo m j_ n , we again have 0 ~ b and then a 
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increases further, up to the angle w r/lax = w min + n * The value 
of cu nin increases with f/5 • For = 0, <o m j_ n = 0, hence 
the value of a is smallest on the trailing edge and greatest 
on the leading edge. Conversely, for — - 0, w m in = (tt/e) - A 
and is therefore situated in the middle of the lower surface- 
In general, with increasing f/6, the minimum value of a 
moves from the trailing edge to the middle of the lower sur- 
face; the point where o ~ again from the trailing edge to 
the lea-ding edge, while the maximum value of o moves simul- 
taneously from the leading edge to the middle of the upper 
surface* 

In order to get an idea of the course of the velocity, we 
must now estimate the mutual effect of the factors F and o 2 . 

I will proceed with this discussion in close connection with 
the diagrams, which I must first explain- Their arrangement 
is the same as for the diagrams in Eiffel's "Resistance dc 
Pair." Each figure has, at the bottom, an accurate outline 
of the wing section. Vertically above each point of the wing, 
there is plotted from a zero line on the vertical the ratio 
q 2 /V 2 , the upper curve corresponding to the upper side and 
the lower curve to the lower side of the wing section. The 
dashed line shows the unit distance from the zero line. The 
area enclosed by the q a /V 2 curve gives, when multiplied by 
— V 2 , the lift of a unit width of the wing. .The chosen an- 

<dg 

gle of attack is 6°(6 = C.l), the air flow being horizontal. 
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I am dividing the discussion into several paragraphs. 

1. The strong suction on the greater portion of the upper 
surface is common to all the figures. This is indeed the chief 
source of the lift, while the pressure on the lower surface 
contributes only a small increment. This can be easily veri- 
fied from the general laws. In fact, as already stated, F 
diminishes very slowly along the under surface from the tread- 
ing edge almost to the entering point. Hence F differs -but 
little* on most of the lower surface, from the value 

cos (A + (2), which it has on the trailing edge, and therefore 

a 2 

only a little from unity. Since also the factor — — falls 

r l ~ 

Cj 

only slightly below 1, up to the vicinity of the entering 
point, q/V is certainly not much smaller than one and hence 
there is only a slight pressure. 

2. Fig. 2 (f = 0 and j = indeed shows a suction 

effect along a portion of the under side. Since the values 
of f arc here less than unity, such a suction effect can 
only be very small. Its appearance is due to relatively large 
values of o 2 and depends essentially on the ratio f : 6. 
Between the trailing edge and w = 2u> Tiin no suction can otour, 
because a is here smaller than lj 2. Any suction effect can 
therefore be expected for only small values of f : 6, where 
w nin is small. For 5=0 any suction effect is entirely im- 
possible, since then corresponds to the leading edge- 
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The suction effect has also been expo 
Eiffel on the ring 11 en aile d'oiseau" 
all the surfaces tested by him can he 


rimer. tally determined by 
which probably alone of 
coroared with a Joukovr; 


wing section.* 


3. Even on the upper side, the course of the velocity is 


characteristically aff 

ected by t\ 

he ratio 

ly shown by Figs. 2-4. 

f/6 is 

express o< 

in the position of the 

angle m 

£ 

■H 

3 

1! 

r*. 

has its maximum value. 

The fact 

that o 


f : 5. This is clear- 
on the upper surface 
+ tt, for which a 
ontinues to increase 


from the leading edge as far as au ax causes the maximum value 


of 

q 

to move farther from the 

leading edge than 

the maximum 

of 

F 

and the fall in velocity 

to be less ra.pid- 

7/e differ- 


entiate "slightly cambered" wings (f/5<3), in which w rflaX in 
near the leading edge, and "highly cambered" wings (f/5>3). 


in which a> rnaX lies nearer the middle of the uoper surface* 
Fig- 2, with f = 0, is a typical example of a slightly cam- 
bered wing. Here the maximum value of a is situated in the 

leading edge and the values of F and a therefore decrease 

* Eiffel, "Resistance de 1’air," 1911, Table XII; also p. 105 
and "Complement," p. 102 ("aile ITicuoort"). The Eiffel fig- 
ures show that the suction effect increases on the under side 
with decreasing angle of attack. This is in agreement with 
our theory, for the factor F increases, as shown by formula 
6, at cvpry point on the under side with decreasing |3 . An- 
other suction effect, which Eiffel finds on the trailing edge 
of nearly all wings, is doubtless due to the formation of 
vortices. 
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simultaneously, thus producing a very pronounced maximum veloc- 
ity, although the maximum velocity is not important in itself. 
On the other hand, Fig. 4 (f/l “ l/5, b/l ~ l/20) shows, in 
spite of a twice as large maximum velocity, a remarkably slow 
velocity decrease toward the upper side- In fact, in these 
experiments, the maximum of a is situated at about l/3 of the 
upper side and a more rapid velocity decrease accordingly 
first begins behind this point. We note also the small inter- 
mediate maximum on the upper side, which is caused by the in- 
crease of a z in spite of the simultaneous decrease in F. The 
mean between Fig. 2 and Fig. 4 is held by Fig. 3, with 
f/l = l/lO and b/l ~ 1/20. Here does not lie very far 

from the leading edge, about l/6 of the upper side, the in- 
crease in ch 3 vanishes under the decrease in F and along the 
greater portion of the upper side we note a uniform falling 
off in velocity, due to the simultaneous decrease in the fac- 
tors F and o u 

These relations were also found in Eiffel's experiments 
with the wing " en aile d'oiscau." Even the intermediate maxi- 
mum of q on highly cambered wings is found on his figures.* 

Lastly, I wish to call attention to the fact that the 
measurements of Fig. 4 appear to me to be worthy of commenda- 
tion, on -'account of the very uniform stressing of the upper 

side. 

* On Eiffel's figures, it appears that, with decreasing p , 
the maximum velocity moves backward from the leading edge on 
the upper side. This also agrees with the theoretical conclu- 
s ions". " 
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4. Fig. 5 has a very slight rounding (6/l = l/50) at 
f/l = l/5.' Therefore k t = 2,5 ana hence the very high ve- 

o O 

locity maximum on the leading edge. IVe have already seen that 
the high maxima must decrease very rapidly toward the upper 
side- The region of this steep decline corresponds to an angle 
of about the size (3 . During the drop, however, there is in 
the figure a long space of almost constant velocity. This is 
explained, as in paragraph 3, by the fact that the maximum 
value of a is located at about 1/2.5 of the upper side- Only 
behind this point is there again a rapid decline to the trail- 
ing edge. This behavior is generally characteristic for high- 
ly cambered wings of slight rounding and occurs also on Eif- 
fel’s diagrams. 

As regards the production of the diagrams, it may be noted, 
in conclusion, that they were drawn according to the very con- 
venient method of E- Trefftz, s.s set forth in the accompanying 
note. Wherever it appeared necessary, the plotting was veri- 
fied by calculation. 
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